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I. INTRODUCTION 


Heavy baryons combine relativistic light quarks and non-relativistic heavy quarks and may have similarities with heavy 
light mesons and, in the case of doubly charmed baryons, also quarkonia. These particles can be understood using a number of 
theoretical tools, in particular. Heavy Quark Effective Theory (HQET) |[T]|3, non-relativistic QCD (NRQCD) fSHU and potential 
NRQCD (pNRQCD) ||5]|6l. Currently, 19 charmed baryons and 8 bottomed baryons are present in the PDG summary table^jTl . 
Shortly after the discovery of the the first charmed baryon was detected in 1975, the A+, at the BNL ||8l, followed by the 
discovery of E++ in 1976 at ENAL and the first bottomed baryon, the Af, in 1981 at CERN IfTOll . In the last decade, bottomed 
baryons were studied at the Tevatron and more recently at the LHC, whereas charmed baryons were mainly discovered at the 
B-factories. While masses, lifetimes, widths and form factors have been measured, direct spin and parity identification is often 
still missing (they are assigned from quark model considerations). This situation will probably improve with the study of angular 
distributions of particle decays in LHC data thanks to the large statistics. In the future, the Belle-II experiment at SuperKEKB 
and the PANDA experiment at the PAIR facility will study singly charmed baryons. With respect to doubly charmed baryons, the 
prospects are not so promising. Only SELEX published evidence for the (3520) ifTTlfT^ and its isospin partner S++ (3460) 
im, with no further confirmation by any other experiment. Besides the experimental findings, various theoretical approaches 
have been employed. Among these are studies based on quark models iflT^E^ . HQET QCD sum rules ll24^j^ and lattice 
QCD, using the quenched approximation Il30l - l34l and, more recently, including sea quarks ll35H44l . 

While most previous studies have focused on post- and prediction of the charmed baryon spectrum, in this work we investigate 
the light flavor dependence of the singly and doubly charmed states. Observed spectra of mesons and baryons have long been 
understood in terms of flavor symmetry, with, for the example of three light quark flavors (Nf = 3 corresponding to up, down 
and strange), the mesons falling into singlets and octets and the baryons into octets and decuplets. While SU(2) isospin symmetry 
is reasonably well respected in nature, SU(3) flavor symmetry is not. The pattern of symmetry breaking for the latter can be 
derived by expanding about the flavor symmetric limit in the strange-light quark mass difference {6m£ = mg — me). This leads 
to the Gell-Mann-Okubo relations ll45l l46l . which are found to hold within a few per cent. Enlarging the symmetry group to 
SU(4) to include the charm quark provides indications of what charmed mesons and baryons should exist, however, the mass 
spectra may be best explained by treating charm quarks as spectators and considering the SU(3) symmetry breaking pattern as 
for the light hadrons. 

In the past, studying the flavor structure of hadrons on the lattice was mostly restricted to approaching the physical point, 
keeping the strange quark mass approximately constant (and consistent with experiment) while reducing the up/down quark 
mass. Extrapolations to the physical point, if necessary, are guided by chiral perturbation theory, which is of uncertain validity 
in the range where lattice results are generated depending on the observable Il47l . Recently, the QCDSE collaboration adopted 
an alternative strategy ll48l l49l where, starting from the Nf = 3 theory, one approaches the physical point keeping the average 
quark mass fixed. This enables one to derive the quark mass dependence of physical quantities a la Gell-Mann-Okubo. The 
extrapolation to the physical point is limited by the order of the expansion required relative to the number and precision of the 
data points available to fix the corresponding coefficients. The incorporation of an additional valence quark, for example, the 
charm, into this framework is described in ll49l . In this work, we study the flavor dependence of singly and doubly charmed 
baryons at a single lattice spacing for both positive and negative parity states treating the charm quark as a spectator, with the 
view to performing a larger scale analysis on CLS ensembles iBOll including a continuum extrapolation. Eor the associated 
Gell-Mann-Okubo relations to be applicable, the charmed baryons must fall into the expected SU(3) multiplets. With only a few 
exceptions we find this to be the case. The SU(4) representations naturally suggest interpolators for creating (and destroying) 
the baryon states. We also compare these interpolators with another basis derived from HQET. Similar results for both bases 
are obtained for positive parity states, while differences are found for negative parity states. Eor the latter, some channels are 
obscured by the presence of two-particle scattering states of the same quantum numbers. We carefully study this as well. 

The article is organized as follows. A description of our computational setup is given in Section]^ This includes a discussion 
of the leading O(a^) discretization effects, which could be significant as the simulations were performed on ensembles at a single 
lattice spacing a ^ 0.075 fm for which the charm quark mass in lattice units is around 0.4. We attempt to quantify these effects 
by comparing results for the lower lying charmonium and D/Dg states to experiment. In Section III we motivate our basis of 
interpolators for the correlation functions, followed by a discussion of the methodology used to extract the mass spectrum and 
the efficacy of our interpolator basis in Section |IV] The charm quark is partially quenched in this study (it does not appear in 
the sea) and as such must be tuned to reproduce experiment on each ensemble. This procedure is described in Section [V] We 
remark that for transparency and in order to maximize the predictive power of our simulations, we try to make use of as little 
experimental input as possible. Eor instance, we predict the absolute meson and baryon masses using one and the same value of 
the charm quark mass, rather than quoting splittings relative to a reference mass, such as that of the rjc meson. Comparing two 
volumes, we also investigate finite volume effects in SectionP 


VIA 


point of the charmed baryon spectrum is given in Sections 


VIB 


The main result of the paper, the extrapolation to the physical 
and VIC and for the mass splittings in Section [VID| before 
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comparing with other recent lattice determinations in Section VIE We finish with some concluding remarks in Section VII 


Additional details are provided in Appendices [A] (meson effective masses), [^(finite volume effects), [C| and [P] (derivation of the 
Gell-Mann-Okubo expressions) and (fit ranges and extrapolations). 


II. SIMULATION DETAILS 


We have employed SLiNC (Stout Link Non-perturbative Clover) ifSTl Nf =2 + 1 gauge configurations, generated by the 
QCDSF collaboration BSl |4^ . The gluonic action is tree level Symanzik improved and the fermionic action has a single 
level of stout smearing in the hopping terms and unsmeared links in the clover term. The clover coefficient was determined 
non-perturbatively. The quark masses were chosen by finding the SU(3) symmetric point where the average octet pion mass, 
= -y/(M^ + 2M^)/3, coincides with experiment. and Mk correspond to the pion and kaon masses, respectively. 
The strange and light sea quark masses are then varied so as to approach the physical point keeping the singlet quark mass, 
TO = |(2 to^ + TOs) fixed up to 0(a) corrections. The ensembles used in this work include several pion masses and two 
volumes for a single lattice spacing, see Table Note that the physical SU(3) symmetric value = Mk « 411 MeV was 
somewhat missed. Below, we discuss our strategy for correcting for the “wrong” trajectory in the mass plane and also quantify 
the size of discretization errors in spectral quantities. In order to reduce auto-correlations, a single measurement was performed 
per configuration where the position of the source was randomly chosen and consecutive configurations are separated by two 
trajectories. 


KE 



Kc2 

L/a X T/a 

(MeV) Mk (MeV) LM^ 

^meas 

(fm) 

1 

O 

flis-ip (fm) 

0.12090 

0.12090 

0.11065 

0.1116 

24 X 48 

471 

471 

4.3 

2747 

0.0756(10) 0.076(2) 


0.12100 

0.12070 

0.11065 

0.1116 

24 X 48 

394 

493 

3.6 

1018 


0.076(2) 


0.12104 

0.12062 

0.11065 

0.1116 

24 X 48 

364 

507 

3.3 

934 


0.076(2) 


0.12090 

0.12090 

0.1110 

0.1116 

32 X 64 

461 

461 

5.6 

875 

0.0748(9) 


0.072(4) 

0.12104 

0.12062 

0.1110 

0.1116 

32 X 64 

355 

499 

4.3 

989 

0.0742(8) 


0.068(4) 

0.121145 

0.120413 

0.1110 

0.1116 

32 X 64 

259 

530 

3.2 

885 

0.0746(9) 


0.075(4) 


TABLE I. Details of the ensembles used, acs and k,ci/c 2 correspond to the light (u/d), strange and charm quark mass parameters, 
respectively, and the gauge coupling is the same in all cases, /3 = 10/g^ = 5.5. Nmeas indicates the number of measurements used for the 
analysis. The last three columns indicate the values of the lattice spacing for the different ensembles, determined through wq, ro and the IS 
- IP charmonium splitting. Note that the lattice spacing determinations were performed using different numbers of configurations than those 
used in the spectroscopy analysis. 
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FIG. 1. The lattice spacing determined via wo, ro and the IP-IS charmonium splitting for each ensemble. 

The valence charm quark is treated relativistically using the SLiNC action. We employed two different values of the mass 
parameter for each set of configurations, given in Table [I] These were chosen so that we could interpolate to the physical charm 
quark mass determined by comparing the spin-averaged IS charmonium mass to experiment (see Section|V]l. In order to convert 
dimensionless lattice results to physical units the lattice spacing must be determined. We considered three quantities to set 
the scale; the spin-averaged IP-IS charmonium splitting, the Sommer scale, rg ll52l (where we used rg = 0.50(1) fm) and the 
Wilson flow observable wg ll53l . which is related to tg E). Fig.[^shows the results for each quantity, where available. The values 
obtained are reasonably consistent across the ensembles and also from the different observables. We take a = 0.075 fm without 
quoting an error|^ since the systematics due to the lack of a continuum limit extrapolation will induce a larger uncertainty. 


^ Note that we take the continuum limit value of wq at physical quark masses. Therefore the precise lattice spacing is ambiguous up to a few percent. 
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Our value of a is about 10% smaller than the value given in Refs. Il48l l49l (a ^ 0.083 fm) which was set from the flavor 
singlet baryon mass QCDSF obtained. However, it is consistent with a later determination by QCDSF ll55ll (a ^ 0.073(2) fm). 
Note that the smaller lattice spacing means our values for the pion masses given in Table [I] differ from those of Refs. Il48ll49l . 
In addition, the average octet pion mass, X^r, is larger than the experimental value by approximately 50 MeV (V = 32^ x 64). 
Recall that m is kept fixed in our simulations and = XT^irn) is approximately constant as Smi = rus — mi is varied from the 
flavor symmetric point {5mi = 0) to the physical value where ~ « 0.225 

GeV^. This means that an extrapolation to the physical pion mass will result in an unphysically heavy kaon, as illustrated in 
Fig. 1^ (point C). SU(3) mass multiplets are extrapolated via Taylor expansions in 5mi starting from the flavor symmetric point 
such that an individual hadron mass has the dependence M = M^m) + c6mi + 0{6m'j). Note that the linear coefficient of 
this expansion (c) does not depend on m and, up to quadratic corrections, Mo(to) « Xmuitipiet- The latter quantity denotes 
the flavor average for the given multiplet. In order to make contact with the physical theory we extrapolate masses within each 
SU(3) multiplet to (point B in Fig.j^ and then shift these masses by an estimate of (i.e. moving 

from B to A in Fig.[^. This procedure is described in more detail in Section 


VIB 



Mi {GeVf 


FIG. 2. The pion and kaon masses squared in physical units, using a = 0.075 fm, for our ensembles (blue circles) compared to the physical 
trajectory. The point A indicates the physical pion and kaon masses, while B indicates the point on the simulation trajectory corresponding to 
the physical value for Smi = rUs — me and C shows the (unphysically heavy) kaon mass corresponding to physical pion mass (but unphysical 
Sme). 


Particles 

Operators 

r/e, D, Ds 

517592 

J/ii,D*,D*s 

5i7i52 

XcoYP), Dq, Dso 

5152 

Di, Dsi 

51717552 , 5i57'fc7i7fc52 

Xcl(lP) 

51717552 

heilP) 

qitijkXjXkq2 


TABLE II. The meson interpolators employed to calculate the charmonium, D and Ds spectra. 


For the SLiNC action the leading discretization effects are of O (a^). The combination O ((mco)^) may be significant for 
some quantities given that am^ 0.4 in our simulation. In order to gauge the size of the corresponding systematic uncertainty 
we have computed the low-lying charmonium, Dg and D spectra for the V = 32^ x 64 ensembles using the interpolators given 
in Table Correlators constructed from these interpolators with quarks smeared over different spatial extents were combined 
with the variational method (see Section |IV]i to extract the ground state. The extrapolation to the physical point is discussed in 
Sections [VI B and VIC the results are presented in Fig.and Table |nl| 

Overall the radial and orbital excitations are reasonably well reproduced. This is to be expected since the typical energy 
scale for the IP-IS and 2S-1S splittings in heavy-light systems is A ~ 0.5 GeV, which is much smaller than the inverse lattice 
spacing. Similarly, for the charmonium ground state the corresponding energy scale is rngV^ ^ 0.5 GeV for ~ 0.4 ll5^ . 
We find for charmonium that the r]c{2S), '0(2S), Xci and he states lie 50 — 60 MeV below experiment corresponding to 2 — 3.7 
standard deviations. However, as discussed below, the 2S and IP fine structure splittings are reproduced, albeit with large errors. 
It is likely that other systematics, in particular, finite volume effects are important for these radial and orbital excitations. These 
systematics are not investigated for mesons in this paper, however, finite volume effects for the charmed baryons are discussed 
in Section IVI Al 
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FIG. 3. The low lying D (top), Da (middle) and charmonium (bottom) spectra at the physical point for the V = 32^ x 64 ensembles. See 
Sections[VI B|and|VI C|for details on how the results at the physical point are obtained. 


Channel 


M (GeV) 

Channel 


M (GeV) 

Vc 

0 -+ 

2.9929(12) 

7141968) 

0 " 

1.9824(85) 

rye (2s) 

0 ”+ 

3.5778(296) 

7i:(2112) 

1 " 

2.1054(95) 

J/ip 

I" 

3.0944(19) 

714(2317) 

0 + 

2.3490(191) 

iP{2s) 

1 ” 

3.6294(246) 

llsi(2460) 

1 + 

2.4415(271) 

XcO 

0 ++ 

3.3914(228) 

71si(2536) 

1 + 

2.5092(431) 

Xcl 

1 ++ 

3.4570(192) 

D 

0 " 

1.8778(106) 

he 

1 +- 

3.4697(150) 

D* 

1 " 

2.0041(142) 




Hi 

1 + 

2.4147(610) 


TABLE III. The low lying open and hidden charm meson spectra at the physical point for the V = 32® x 64 ensembles. 

For the 0+ and 1+ heavy-light mesons that are above or close to strong decay thresholds one needs to consider the relevant 
scattering states. We are able to resolve two closely lying states for the Dsi by using two interpolators (see Table [IJ in addition 
to multiple smearings in the variational method. However, a proper finite volume analysis would be required to identify the true 
nature of the higher lying state as the D*K threshold lies between the D^i (2460) and the Dsi (2546). For the Di the same analy¬ 
sis showed the lowest eigenvalue of the variational method to be clearly consistent with D*tt, and the next level to be compatible 
with experiment (the latter shown in Fig.[^, higher eigenvalues were much larger in mass. As for the Dsi, from heavy quark 
symmetry one expects two states close together Il57l - l59l and in experiment there is the 11°(2420) (width 27.4(2.5) MeV ITJ) and 
the 11°(2430) (width 384)*)}°° MeV Q). However, a larger basis, including interpolators with derivatives, would be required to 
resolve the additional level. For the 0+, only one state is expected and our results from a single scalar operator are compatible 
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with experiment for the £)*g within 2cr. Exchanging the strange quark for an up/down quark, the lowest eigenvalue of the Dq is 
consistent with Dtt. Unfortunately, we were not able to reliably extract the second eigenvalue for this channel and so we do not 
include a value in Fig.[^ A more extensive analysis of open and hidden charm mesons with a larger set of interpolators, including 
also J = 2 states, will be presented in a forthcoming publication. We note that studies of open charmed meson channels near 
thresholds including four quark interpolators have already been performed, see Refs. ll^lMI . 

Fine structure splittings are dominated by higher energy scales (O(mcVc) and 0(mc) in charmonium and the D/Dg systems, 
respectively) and are therefore more sensitive to discretization effects. We find Mj/^ — — 100.1(1.5) MeV, Mjj* = 

122.0(2.6) MeV and Mjy — Mu = 126.8(10.1) MeV. These splittings are approximately 13 MeV, 22 MeV and 15 MeV below 
experiment, respectively. Assuming lattice spacing effects are the main cause for the discrepancies, we estimate 10 — 20 MeV 
as the likely size of this systematic in the charmed baryon spectra (in particular for spin splittings]^ For most observables 
the total error of our hnal results after extrapolation to the physical point is of a comparable or larger size. For the radial and 
orbital excitations in charmonium, the wave function will be broader, suggesting smaller lattice spacing effects. We find the 2S 
hyperhne splitting, ip — ric = 48(20) MeV and the IP splitting, Xci — he = 11(19) MeV, compared to 46.7 MeV and 14.8 MeV, 
respectively, in experiment. 


III. CHARMED BARYONS AND INTERPOLATORS 


The simplest way to see which charmed baryons are likely to exist is to consider the irreducible representations of the tensor 
product of three SU(4) fundamental representations: 


4 (g) 4 (g) 4 = 20s © 20 m © 20m © 4a. 


( 1 ) 


This flavor symmetry is not respected in nature, however, the number of baryons and their flavor quantum numbers may be 
reproduced. Fig. displays the totally symmetric 20s-plet, the mixed symmetry 20M-plet and the total anti-symmetric anti¬ 
quadruplet, 4A- For each state within a multiplet the total spin and parity is the same and resonances may also be expected to 
fall into this pattern. The ground state 20s and 20M-plets have positive parity (consistent with the light baryon members of 
the multiplets) and J = \ and J = 4, respectively, while the A a ground states, in the non-relativistic limit, require non-zero 

orbital angular momentum and have negative parity (J^ = 4 ). Orbital angular momentum is also needed for the negative 
parity counterparts to 20s and 20m- The lowest mass members of the SU(4) multiplets are the familiar decuplet, octet and 
singlet of SU(3) symmetry. The singly charmed baryons can form SU(3) sextets and anti-triplets (3 (g) 3 = 6s + 3a), where the 
restriction of a totally anti-symmetric wavefunction means that there is no anti-triplet within the 20s-plet. The doubly charmed 
baryons are all in SU(3) triplets. 

In nature, where SU(3) flavor symmetry is broken, the physical states may not fall exactly into the expected multiplets, 
for example, the Sc (S^,) may contain a small contribution from the sextet (anti-triplet). Furthermore, the above picture does 
not include possible non-quark-model states. However, the observed charmed baryon spectra, with either measured or 
assignment from potential model predictions (for a review see, for example. Ref. ||62]) seem to reproduce the expected pattern: 
replacing a strange quark by a light quark, splittings are similar within a given SU(3) multiplet, consistent with a constituent 
quark model picture. We remark that since we are simulating QCD in the isospin limit (m„ = nrid and omitting QED effects) 
then all states within an isospin multiplet are degenerate. 


^ We note that there is an additional uncertainty arising from the scale setting of a few percent in the splittings. 
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FIG. 4. The SU(4) irreducible representations for baryons: (left) the totally flavor symmetric 20s-plet, (middle) the mixed symmetry 20M-plet 
and (right) the totally anti-symmetric anti-quadruplet. 


The internal structure of baryons containing light {q G {£, s}) and heavy quarks (Q) can be addressed in terms of HQET and, 
for the QQq case, pNRQCD, starting from the static limit. At finite heavy quark masses also NRQCD is valid for doubly heavy 
systems. In the case of Qqq baryons, the heavy quark Q provides a color source for the two light quarks. In the mq —> oo limit, 
the light quarks have a definite total angular momentum and a total spin = 0 or 1 , corresponding to a flavor anti-symmetric 

or symmetric structure, respectively. In this limit, the spin splittings between baryons with Sd = 1 and = 0 vanish since they 
_ 2 _ 

are of 0(A /mq), where A is the energy scale of the light degrees of freedom. 

For doubly heavy baryons, in what we call the HQET picture, the two heavy quarks form a diquark of small spatial extension 
that interacts with the light quark in analogy to heavy-light mesons (if the QQ diquark is in a color anti-triplet), shown in Fig. 
HQET corresponds to pNRQCD in the limit of the distance between the two heavy quarks, r —> 0. Assuming such a QQ 
diquark, to leading non-trivial order in Ijmq, one can show that QQq baryon spin-splittings, Mqqq{J = |) — Mqqg{J = i). 


are 3/4 times the corresponding Qq fine structure splitting, Mq^{J = 1) — Mq^{J = 0) 

Alternatively, the doubly heavy baryons could be comprised of a heavy and a light quark in a (color-anti-triplet) diquark 
and, together with the remaining heavy quark, one has a charmonium-like system, also shown in Fig. In this case the level 
splittings can be understood in terms of pNRQCD and NRQCD, for example, the QQq spin splittings can be related to the 
charmonium fine structure splitting. However, this is not so straightforward as for the HQET picture since the light quark within 
the Qq diquark cannot be considered as spatially localized. It is possible that the HQET picture works best for the lower lying 
states (where r <C A ), while a charmonium-like picture is applicable for higher excitations with r > A . pNRQCD includes 
both possibilities. 

The expected internal structure of a particle informs the choice of lattice interpolator employed, since we want to have a good 
overlap with the physical state. Any interpretation that is valid in the heavy quark limit may not work particularly well for 
charm quarks and will be more applicable to baryons invol ving bottom quarks. With this in mind, we implemented two sets of 
interpolators, those based on HQET for J = ^ and | (Table llVl where the diquark is formed from two light quarks or two heavy 
quarks in the case of singly and doubly charmed baryons, respectively) and also those arising from SU(4) symmetry (Table [V|) 
for J = i. A comparison of the efficacy of these interpolators is made in Section 


IV 




FIG. 5. Internal structure of a doubly heavy baryon: HQET picture (left hand side), Quarkonium-like (right hand side). The blue circles 

represent the heavy quarks, Q, aiid r is the average Q — ©.separation. Jn the HQET picture r <C A , while in the charmoniqmpictuire r > A,. 
The mterpolarors given in the tables do not nave dennite parity and one tfeeds to use the projection oj 3 erator ± 74 ) 

to obtain positive or negative parity states. Since, in the non-relativistic limit, negative parity requires non-zero orbital angular 

momentum, interpolators including derivatives may improve the overlap with the physical states. Further exploration of the 

best basis will be performed in a future study. An additional projection is required for the HQET interpolators of type and 
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HQET singly charmed baryon interpolators 


S I Sd 

{qq)Q 0 

J = h J = l 

0 0 0 

{££')c Os-y = ea6c(r^C75f'*')c;( 

A, 

0 1 1 

{££)c 

Ec E* 

-1 1 0 

{£s)c Osj = eabc(£°''^C'y5s'’)Cy 

“C 

-11 1 

(£s)c = £ate(r^C' 7 ^s *’)4 

“C ' —'C 

-2 0 1 

(ss)c = eate(s“'^C' 7 ^s*’)c^ 

He H* 


HQET doubly charmed baryon interpolators 


0 1 1 

{cc)£ Oy^y = ee 6 c(c“^C 7 ^,C*’)f= 

dec “CC 

-1 0 1 

(cc)s Oy^y = eaf,c(c“^C7Mc‘')s" 

0 0 * 

iicc 


TABLE IV. Quantum numbers of singly and doubly charmed baryons and interpolators, Oa-^ in the HQET picture, where 7 is the spin index. 
i and I' stand for up and down quarks, c for charm and s for strange. S, I and Sd are strangeness, isospin, and diquark total spin quantum 
numbers, respectively. 


SU (4) singly charmed baryon interpolators 

SU (4)-plet S I 

0 


20m 0 0 

Ob-y = ^e“*'"|2(r^C'75f2'’)c} -b C-id'‘)ty - C'^d)£y] 

Ac 

-1 1 

Osy = ^e“‘"{2(s“^C75f|)c;j + - (c“^C75s'’)f" } 

“c 

0 1 

Oby = e“*'"(c“^C'75f'’)f^ 

Ee 

-1 1 

Os-y = + (r^C775c'’)s^} 

'—'C 

-2 0 

Oby = e“*'"(c“^C'75s'’)s^ 

He 

SU (4) doubly charmed baryon interpolators 

20m 0 1 

Os-y = C 75 c'’)c^ 

“cc 

-1 0 

Osy = e“*’"(s“^C75c'’)c^ 

^cc 


O'^ (see Talftl^^ys AesBTdii^^jBottJMiate^QlsifctJ feo^&Htfparyeiiinetil/hfetUitscefaralsialictniBPbcei^sfeffltinglltLtiStag projection 
operators (at zero momentum): 


= 5^ - L j e {1,2,3}, 

( 2 ) 

We note that the interpolators for the and states will, in principle, have an overlap with both physical states since they 
are not distinguished by a conserved quantum number. In the HQET case, the two interpolators differ in terms of the total spin 
of the diquark {sd) which is not conserved at a hnite heavy quark mass and for the SU(4) case the interpolators are in different 
SU(3) flavor multiplets, however, the flavor symmetry is, of course, broken. We have not taken possible mixing between the 
states created by the Sc and S} interpolators into account. This would require calculating the corresponding cross-correlation 
functions. However, such a study was performed in Ref. Il43l with high statistics and no significant mixing was found. This 
suggests that the approximate SU(3) symmetry holds sufficiently well to suppress the mixing. 


IV. VARIATIONAL METHOD AND FITTING PROCEDURE 

We compute the mass spectrum in the conventional way by calculating two-point correlation functions created from sets of 
baryonic interpolators with spin index 7 . These correlation functions contain contributions from all states with the same 
quantum numbers given by the interpolators 

[C{t)],, = = ^(0|0.,-,(0)|n)(n|0,,^(0)|0)e-^"‘, (3) 

n 

where is the parity projection operator. We employ the variational method llMll65]l in order to reliably extract the ground 
state and first excited state. For each interpolator given in Section [111] we generate a 3 x 3 matrix of correlators by varying the 
smearing applied to the quark fields. Extracting the eigenvalues, A“(f, fg), of the generalized eigenvalue problem at large times, 
the lowest two states are cleanly separated: 




( 4 ) 
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where AMa'a = Ma' — Ma, and a' > a. The time-slice to can be varied to reduce the excited state contaminatiorj^to A“(i, to) 
and the eigenvectors to)- We use Wuppertal smearing ll^ l68l with APE smoothed links ll^ iTOl where the number of 
Wuppertal iterations is varied to obtain one smearing combination which leads to a spatially extended interpolator with a good 
overlap with the ground state and two other combinations which have significant, but different, overlap with excited states. For 
the singly charmed baryons, the overlap of the interpolator with the physical states was most sensitive to the light quark smearing. 
The number of Wuppertal iterations for the heavy quark was fixed to n = 150 and for the light quark n G {5,25,150} was 
realized. For the doubly charmed baryons the situation is reversed and n = 150 was used for the light quark and n S {5, 25,150} 
for the heavy quarks. 



FIG. 6. Effective masses for the first (EVl) and second (EV2) ei^nvalues of Ecc for positive and negative parity obtained using different 
values for to when solving the generalized eigenvalue problem Eq. for SU(4) interpolators on the symmetric ensemble with V = 32® x 64. 


The eigenvalues A“^®’^(t,fo) were fitted (separately) with single exponential^ for fixed tg > 2a in the range ti to tf, 
taking the correlations between time-slices into account. The third eigenvalue was discarded as M 3 cannot be cleanly separated 
from higher excited states. The final fit ranges chosen, compiled in Table |XI| (Appendix [^, have reduced correlated values 
X^/dof < 2 and in most cases x^/dof ~ 1. For these fit ranges the masses extracted were stable within errors as ti was further 
increased. No significant dependence on tg was found (see, for example. Fig. and we take tg = 2a. The statistical errors 
were evaluated using the jackknife method combined with binning. Measurements were performed on every other trajectory for 
each ensemble and the errors were stable for nbin >2 — 4. We made the conservative choice of ribin = 4, which is consistent 
with Tibin > dTint, where Tint is the integrated auto-correlation time. The latter was estimated via the F-method 117111721 to be 
between 0.5 and 0.7, depending on the state. 

Effective masses. 


M^{t) = 



/ X°^{t-a,to) \ 
\X°{t + a,to)) 


(5) 


and fits for a representative sample of states are shown in Fig. [7]for the HQET interpolators on the symmetric ensembles with 
V — 32^ X 64 corresponding to M^r = Mk = 461 MeV. For the positive parity states we are able to extract a reasonable signal 
for the first two eigenvalues for both J = ^ and |. For the negative parity states, which are statistically noisier, clear ground 
state signals were obtained for J = \ and |, however, the first excited state was only reliably extracted for the doubly charmed 
baryons. 

For the spin-1/2 particles we can compare the HQET and SU(4) interpolators. Effective masses for a state from each multiplet 
are shown in Fig.|^ The positive parity particles display a consistent picture: the SU(4) interpolators have a slightly better overlap 
with the desired states for the sextet baryons (Sc, S}, Hc)j are only marginally better for the anti-triplet baryons (Ac, Sc) and 
no difference is observed for the doubly charmed triplet (Sec, Hcc). The latter is due to the corresponding correlation functions 
only differing by terms which are suppressed in the non-relativistic limit, when expressed through upper and lower components 
of the quark spinors. This is not the case for singly charmed baryons nor for negative parity. Indeed, a more striking pattern 
emerges for negative parity. For the sextets, the SU(4) interpolators have a much better overlap with the ground state while for 
the triplet and (even more markedly) the anti-triplet the HQET interpolators are clearly better. For the first excitation the SU(4) 
interpolators provide a better signal, although there is no clear single state dominance for the singly charmed baryons. 


'* For t > to > t/2, only a' > N will contribute where N is the rank of the basis 1661 . a limit that we do not consider here. 
^ We assume the baryon masses to be heavy enough for the backward propagating paificle to have a negligible influence. 
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FIG. 7. Effective masses for the first and second eigenvalues for HQET interpolators on the symmetric ensemble with V = 32^ x 64. The 
filled regions indicate the fit ranges chosen and the fit results, including the statistical errors. 


These comparisons are made with the same smearings applied to both sets of interpolators. Changing or optimizing the 
smearing in the individual cases may change these conclusions. The SU(4) and HQET interpolators in each channel both 
belong to the same (lattice) hypercubic representation; either the dimension two Gi representation for the interpolators denoted 
as spin-1/2 in Tables IV and IV] which in the continuum gives J = i, |,... or the dimension four H representation for the 
spin-3/2 interpolators corresponding to continuum J = In the limit of very large statistics the SU(4) and HQET 

interpolators should give the same energy for each state. Eor positive parity consistency is seen, while for negative parity we use 
the interpolator with the best overlap for our setup. Similar behavior is seen throughout for the smaller volume, V = 24^ x 48. 

Note that care has to be taken in order to avoid the misidentification of multi-particle states, for example, a positive parity 
baryon plus a pion as a negative parity baryon. We have studied this systematically for the ground state negative parity channels. 
All but one of our = | ground states are smaller in mass than two-particle states of the same quantum numbers, consisting 


71 + 1 + 

either of a | baryon or a ^ 


None of the ^ 


states contain D or 


2 baryon plus a pseudoscalar meson in a S-wave or a P-wave, respectively. 

mesons - all combinations of these mesons with baryons of the same combined isospin, strange 
and charmness are heavier in mass. There are, however, quite a few S-wave decay channels with thresholds close to our mass 
estimates that require careful study, in particular 


Ac —>■ Ec + tt , 

Sf; —y A “h TT, 

Yjc —y Sc + TT , 

'zr* V 11'* 1 

^ 4- TT , 

^ 4- K , 

^ H" at , 

: 

^c “c “1" ’’’ J 

“c ■'^c + at , 


[I] 

4 

[I] 

o 

-F 

4 

[I] 

4 

[I] 

-F 

4 

Oc + A', 

ilc ^ K J 

“CC ^ “CC “1” ^ ? 

^cc ^cc + at , 


where the baryons on the left hand sides are implied to have negative parity and those on the right hand sides positive parity. 
Neither have we made the isospin nor charges explicit. Comparing all these channels at our three different sea quark mass 
combinations leads to the following conclusions. The ground state negative parity Qc state, obtained both from the SU(4) and 
HQET inspired interpolators, is always degenerate with the sum of the 5^ and a kaon. We show the corresponding effective 
masses at = 461 MeV and at = 259 MeV in Eig. Likewise, we identify our negative parity signal as a Ac 

plus kaon scattering state. 

Particularly interesting is the behavior of the negative parity Ac, also depicted in Pig. |^for two ensembles; while effective 
masses obtained from the SU(4) interpolator are degenerate with the combined mass of a Ec and a pion, the effective masses 
from the HQET interpolator are systematically lower, at least at our two heavier pion mass points. This suggests the HQET 
interpolator to significantly overlap with a physical state lower in mass than the scattering state while the SU(4) interpolator in 
this case strongly couples to the close-by two-particle state. All the other negative parity singly and doubly charmed baryons 
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FIG. 8. Comparison of the effective masses for SU(4) and HQET interpolators employed on the ensemble at the symmetric point with 
V = 32^ X 64. The ground (first excited) state is indicated by filled (open) symbols. 


seem to be relatively stable under strong decays or, in the few cases where their masses are higher than those of the potential 
decay products, at least the interpolators we employ are insensitive to the presence of these scattering states. A similar analysis 
for the first excitation (for both parities) is challenging due to the number of relevant two- and even three-particle states. 


The lower lying meson spectrum in the charmonium, D and Dg sectors was also determined in order to estimate the size 
of discretization errors (as discussed in Section 0 and to enable the tuning of the charm quark mass parameter (see the next 
Section). Except for the 1+ D/Dg channel, we employed 3x3 correlation matrices and applied the variational method, as 
described above, extracting the ground state mass and in some cases the excited state for each meson. Spatially extended 
interpolators where constructed using n S {5, 25,150} iterations of Wuppertal smearing for the charm quark, the light quarks 
were not smearecj^ Effective masses and fit results for a sub-set of states are disp layed in Eig. in Appendix]^ for the 
symmetric larger volume ensemble. Eitting ranges for all channels are given in Table |5Qi| of Appendix]^ For the I'^hannel, 
two interpolators were implemented (see Table [Till and a 6 x 6 correlation matrix was computed. This enabled us to resolve pairs 
of closely lying states, as illustrated in Fig. [^(Appendix [ a]i. In terms of multi-particle states, the 0+ and 1+ D/Dg mesons are 
above/close to the Dtt/DK and D*Tr/D*K experimental thresholds, respectively. We discuss our results for these channels, 
also shown in Fig.[^ in Section VI B| 


® Note that for meson interpolators without derivatives the smearing can be “transformed” from one quark to the other. 
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FIG. 9. Comparison of negative parity effective masses with the relevant strong decay thresholds for (top) Q.c and (bottom) on (left) the 
symmetric ensemble and (right) the most asy mmetric ensemble for V — 32^ x 64 . 

f^s L / Qj X T j CL 6^charm 

0.12090 0.12090 32 x 64 0.1114801(67) 

0.12104 0.12062 32 x 64 0.1114869(43) 

0.121145 0.120413 32 x 64 0.1114908(48) 


TABLE VI. Final values for Kcharm for each ensemble determined using the IS charmonium mass with full statistics. 


V. TUNING OF ^charm 

The charm quark mass parameter, Kcharm, was tuned by requiring that the IS spin-averaged charmonium mass, Mis = 
jM^^ -F is equal to the experimental value. For this quantity the dependence on the light sea quarks is sub-leading 

and furthermore the associated discretization errors are likely to be reduced since contributions from fine structure interactions 
are removed. The tuning was performed in two stages on each 32^ x 64 ensemble. Firstly, a coarse tuning involving 200 
configurations separated by 2 trajectories determined a preliminary value for Kcharm- Two values of the mass parameter, Kd and 
Kc 2 (see Table 1^, were chosen to bracket Kcharm, spaced closely enough that a linear interpolation in 1/k to the physical point 
would be sufficient. In the second stage the meson and baryon spectrum was determined with full statistics for Kci/c 2 and the 
final value of Kcharm determined from Mis- This mass is displayed in Fig. [TOlfor Kci/c 2 and the results for Kcharm are listed 
in Table [Vl] The meson and baryon masses were then interpolated to the physical charm quark mass. Fig, shows there is 
little dependence of Mis at fixed Kcilc 2 on the sea quark masses. This translates into very similar values for Kcharm across the 
different ensembles. The same values were employed for the smaller, 24^ x 48, ensemble where the light/strange quark masses 
are matched. For the intermediate small volume, an interpolation was performed. 


VI. RESULTS 
A. Finite size effects 

With two volumes available with spatial extents, L = 24a and L = 32a, corresponding to 1.8 fm and 2.4 fm, respectively, 
a study of finite volume effects of the charmed baryon spectrum can be attempted. The same quark mass parameters have 
been used for two sets of ensembles, at the symmetric point (Kgym = 0.12090) and the asymmetric ensemble with the light- 
strange combination (Kasym = 0.12104,0.12062), see Table|^ Differences between the masses extracted from the two volumes. 
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FIG. 10. The charmonium Mis mass determined with Kci/c. 2 - The horizontal line indicates the experimental value. 


AM = M 24 — M 32 , can be computed directly; the results are given in Tables VIII and in Appendix [B] The mass splittings 


are only signihcantly non-zero (where we take 2.5 (t to indicate significance) at the symmetric point since we have much higher 
statistics in this case for the smaller volume. However, similar or larger magnitudes for the difference on the asymmetric 
ensemble (although possibly less signihcant) indicate a trend. In addition, we look for consistency between the two types of 
interpolators, HQET and SU(4), where available and that any difference decreases if a light quark is replaced by a strange quark 
within a multiplet (for the asymmetric ensemble). 

In general, one expects larger hnite volume effects for singly charmed baryons compared to doubly charmed ones, similarly, 
for excited states compared to ground states. Whether the effects can be observed depends on the size of the statistical errors. 
Our results fit into the expected pattern. For the positive parity multiplets, we hnd masses are reduced by about 10-40 MeV 
when changing the spatial extent from 1.8 fm to 2.4 fm for the ground state sextets and the anti-triplet, while the ground state 
doubly charmed triplets are unchanged. The statistical errors for the excited states of the singly charmed multiplets are large and 
no signihcant hnite size effects are visible apart from a small increase for the anti-triplet. However, in the doubly charmed case, 
with much smaller errors, clear differences of around 40-87 MeV emerge, again decreasing the mass with increasing volume. 
The negative parity states are less affected with only the ground state anti-triplet showing a notable decrease between L = 24a 
and 32a. The J — ^ sextet and J — ^ triplet have small positive shifts. Note that the J = | negative parity sextet is not 


included in Table IX as we were not able to extract a reliable signal in these channels for the smaller volume. The hnite volume 


effects for a subset of states are shown in Fig. 11 


An estimate of the size of the remaining systematic in the L = 32a results can be made by assuming the hnite volume 
dependence of the masses is of the form Ml = Moo + cF{L), where Moo is the inhnite volume mass and F{L) = is 

the leading order effect in the asymptotic regime and F{L) = 1/L^ otherwise. This leads to 


AMoo = M 32 -Moo = (M 24 - M32)(F(24a)/F(32a) - I)"! = AMAF, 


( 6 ) 


where AF = 0.4 and 0.6 for the symmetric (/vsym) and asymmetric ensembl e (Atas vm) in the asymptotic case, respectively, 
and AF = 0.7 in the 1/L^ regime. Considering the values of AM in Table 


VIII 


this leads to the inhnite volume masses 
at Mt^ — 355 MeV being, for example, approximately 9-22 MeV smaller for the Sc, 11-14 MeV smaller for the Ac and 
34-48 MeV smaller for the excited state Sec and S^c. At this pion mass with L = 32a, LM^ = 4.3, one would expect the 


lower bounds to apply. These estimates should be compared to the hnal values for the masses given in Table VII discussed in 
Section fVI Cl . 

Since the lightest pion mass for the L = 24a ensembles is quite large, M^ = 364 MeV, we do not attempt to extrapolate these 
results to the physical point and they are not included in any further analysis. 


B. Extrapolation to the physical masses 

An important advantage of approaching the physical point using the QCDSF strategy, starting from the Havor symmetric point 
and keeping the average light quark mass hxed, is that the light quark mass dependence of physical quantities is well constrained 
as long as flavor symmetry violations are reasonably small. The dependence of meson and baryon masses on the symmetry 
breaking parameter, 5mi = rus — was derived in Ref. 1491 and extending this method to include partially quenched quarks 
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FIG. 11. Comparison of the V — 24^ x 48 and 32® x 64 results as a function of for a selection of multiplets exhibiting finite volume 

effects, where 1 — oc 5mt + 0((5m|). Top left and top right, the first eigenvalue (EVl) of the positive and negative parity singly 

charmed anti-triplet, respectively. Bottom left and bottom right, the second eigenvalue (EV2) of the positive parity doubly charmed triplet for 
J = I and I, respectively. 


was discussed in Refs. iiiEa. Here, we follow an analogous approach treating the charm quark as a spectator. We classify the 
charmed baryons according to the SU(3) flavor symmetry multiplets for the one or two light quarks present. In the isospin limit, 
we have the following channels 

• sextet: (E,, 5^, H,) for J = ^ and (E*, 5*, H*) for J = |, 

• triplet: (5,,, 11,,) for J = i and (5*„ H*,) for J = |, 

• anti-triplet: (A,, 5,) for J = |. 

The light quark mass dependence of each multiplet can be derived by considering the matrix element {Bii\M\Bfi). The mass 
matrixM = fnl — ^SmiTs = mHo+Sm(,Hs where Tg is the SU(3) generator and is the baryon multiplet in representation 
R. Hg transforms like the eighth component of the octet representation and since m is kept constant, this is the only part of the 
mass matrix which is changing as we approach the physical limit. Performing a perturbative expansion in 6me we have 

{Bii\M\Bn) = + Sm,{Bn\Hs\Bi,) + ^ + ..., ( 7 ) 

where = fn{Bfi\Ho\B]i) is the mass of the multiplet R in the flavor symmetric limit. In the light baryon sector the 

expansion to first order leads to the Gell-Mann-Okubo (GMO) relation ll45ll46ll . The expressions for the linear terms for the 
charmed baryon multiplets are derived in Appendix [P] Here, we comment that the matrix elements are the flavor singlet terms 
arising from the decomposition of the direct products of R 0 Th 0 R, where Th is the representation of H. This means that 
the number of terms appearing at each order is given by the number of times the trivial representation appears in the direct 

product. At the lowest order (R ® 1 0 R) only one term arises, and for the sextet, triplet and anti-triplet that we are considering 

there is also only one term at first order (R ® 8 0 R). Hence, within each multiplet there is only one coefficient at 0{Sme) and 
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one can extrapolate all states within a multiplet (to this order) using just two parameters. However, at second order the trivial 
representation appears three (two) times for the sextet (triplet/anti-triplet) such that there are the same number of new coefficients 
as particles within a multiplet. 

Below we display the expressions for the light quark mass dependence of the charmed baryons masses derived in Ap- 
pendix[^(Eqs. (|D8|l, (D13 1 and (D15 i). We use the GMO relations for the light mesons 


Ml = {M^f - ‘^DSme + 0{5ml), 

(8) 

M]^ = (M® + ^DSmi + 0{6ml), 

(9) 


and + 2M|-) = (M®+ 0{6nnf) to replace Snii with (X^ — M^)/X^ ignoring 0{Sm^) terms. Such higher 

order terms can be absorbed into the coefficients of the quadratic terms below. 


Sextet, R = 6: 


Triplet, R = 3: 


Anti-Triplet, R = 3: 


Me. = M® - 
Mh' = M® + iAi 
Mn. = M® + |Ai 


X^ 

Xl-Ml 

XI 


Mh.. = M® - 


Xl-Ml 


Mo.. = 


0 x2 

A? - Ml 


XI 


Ma^ = - lC^ 


Mh. = M® + \Ci 


Xl-Ml 

XI 

XI - Ml 

XI 


+ A 2 
+ ^3 
+ A 4 


+ i?2 
+ B 3 


+ C 2 
+ C3 



( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


(15) 

(16) 


Note that for the sextet and triplet there are equivalent expressions for the J — ^ states with different coefficients. In 
Refs. 1491 155 I analogous parameterizations are given (but to higher order in their framework^ The above forms can be ap¬ 
plied to extrapolate our results to the physical point if the charmed baryons are regular three quark states where the charm is 
a spectator and the particles fall into the SU(3) flavor multiplets. As for any expansion the coefficients of the flavor symmetry 
breaking terms should be small and the extrapolation will be well constrained whenever only the first few terms are significant. 

In this study we have the additional complication of needing to correct for the fact that the simulation trajectory misses 
the physical point in the rrig — me mass plane (see Fig. [^. The value of to = to^^^® + Am is off by a small amount Afn 
and this means, in the expressions above, the masses Mg also differ from the physical ones by Mg = -f AMq. The 

linear Smg expansion coefficients should be unaffected. Hence, to make contact with the physical point we extrapolate to 
the physical 6me defined as Ml = and then shift the masses within each multiplet by 

AAmuitipiet = where Amuitipiet ^ is the flavor average for each multiplet. 

This procedure should be correct up to 0{AmSml) corrections, see Eq. (j^. For the singly charmed baryons we could choose 


^ Note that their expansion is given in terms of rriq — m, which, unlike = rris — mg, changes when shifting m. 
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FIG. 12. Flavor invariant quantities as a function of M'^ on the 32^ x 64 ensembles. Note that 1 — M^jX^ oc Smt + 0(5m|). (Top left) 
the light vector meson and (top right) the spin averaged IS heavy-light meson. Fits to a constant and a constant plus quadratic term in Smt are 
indicated along with the extrapolated values (blue diamonds) obtained at the point corresponding to physical Sme = The error on the 

extrapolated results incorporates both fits using Eqs. l |17| l and l |18^ . The experimental results are also shown. Similarly for (bottom left) singly 
charmed baryons and (bottom right) doubly charmed baryons for HQET inspired interpolators. For the charmed baryons, the final values are 
obtained by extrapolation to and then shifting by AX^uitipiet, see the text. Two errors are displayed, showing the uncertainty with 

and without the error on AXmuitipiet included in quadrature. 


to determine AXmuitipiet using the flavor averaged J = states and the corresponding experimental results. However, this 
would reduce predictability and the shift would still need to be estimated for the doubly charmed baryons. Instead we look for 
other multiplets that are likely to have similar AXmuitipiet- Assuming that all the singly charmed multiplets are shifted by the 
same amount, we choose the flavor averaged light vector channel, ^{2Mk- + Mp), to estimate AXe « AXg. For the doubly 
charmed triplet we use the spin and flavor averaged heavy-light combination, ^-F -F 2Mjj + Treating the 

charm as a spectator, the D and Dg channels should form triplets and obey modified versions of Eqs. ([^ and ([g. 


Figure 12 shows the light vector and DjDs flavor singlet combinations as functions of As expected for quantities 


which depend on to second order, there is very little change with the light quark mass. Note also that the extrapolation to 
is relatively short. Fitting to a constant (fit 1) and a constant plus quadratic term (fit 2) we obtain extrapolated values, 
Mave ± AMave, using the following formulae 


= 


1 

- I max ( 
2 yiG{l,...n} 


min {Mi-AM,)], 

J 

(17) 

min {Mi-AMi)], 
ie{i,-n} J 

(18) 


where in this case n = 2, Mi ± AM^ is the result from fit f = 1, 2. The difference of Mave to experiment gives the estimate of 
AXmuitipiet- AMave IS taken as the error in the shift. For the D/Dg case the shift is consistent with zero, —3(8) MeV, while for 
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FIG. 13. The shift AX^uitipiet calculated hy comparing different flavor averaged mass combinations with experiment both in the singly 
charmed sector (blue stars) and the light vector meson sector (shaded band). 


the light vector it is 33(17) MeV. 

Figure [T^ also shows the positive parity charmed baryon flavor singlet combinations. Here, the data are extrapolated to 
and the fit error is calculated as above. The appropriate shift, AXmuitipiet is applied, where the associated error is 
added in quadrature. Encouragingly, the final values are in agreement with experiment for the sextet and anti-triplet, consistent 
with AXniuitipiet being small and of a similar magnitude for the different multiplets. This is emphasized in Fig. 13 which 
compares AAmuitipiet determined from ^{2Mk- + Mp) with the shift obtained using the flavor averaged masses for the sextet 
and anti-triplet and the corresponding experimental values for the singly charmed baryons. 

As for the mesons, there is very little dependence o 5mf, for all flavor average masses in Fig. [T^ Ensembles on the correct 
trajectory, close to the physical point will be needed to quantify the size of the leading, i.e. O(Smj), flavor symmetry violations. 
For our data the coefficients of the quadratic terms are zero within 2.5 a. 


C. Baryon and meson masses at the physical point 

Turning to the extrapolation of the individual baryon masses our strategy is as follows. For a given multiple!, we employ the 
Gell-Mann-Okubo relations with and without the quadratic term, denoted GMOl and GM02, respectively. If the xVdof < 
2 for both fits, the final value for the mass at the physical point is computed using Eqs. (0 and ( [T8] | and performing the 
shift, AAmuitipiet. as discussed abov^ The mass values within a multiple! on a given ensemble are highly correlated. These 
correlations have been taken into account in the fits. Overall the GMO expressions describe the data well, however, in a small 
number of cases either the linear (GMOl) fit produced a x^/dof > 2 or both the linear and quadratic (GM02) fits. If the masses 
do not fall into the GMO pattern, we can extrapolate using the functional forms Eqs. ([T0|-([T6]l but allowing the coefficients of 
the linear terms to differ for each member of a multiple! (i.e. only Mq is constrained to be the same). Such phenomenological 
fits are always possible if m is kept fixed in the simulation. Again we considered these fits including and excluding quadratic 
terms, denoted as less-constrained 1 and 2 (LCl and LC2), respectively. Mave and AMave are calculated as above, where for the 



the i anti-triplet are both around 10 — 20%, albeit with large errors. This may indicate they have a more complicated 


* Whenever x^/dof > 1, the error AMi ^2 is inflated by the factor y/x^/dof. 

^ Note that for the LC2 fit there are as many data points as there are unknown coefficients, and thus there are no remaining degrees of freedom. 
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internal structure. We also remark that, as discussed in Section [rv} our ground state spin-1/2 flc and negative parity results 
are likely to be scattering states. In the sextet, only the results for Sc away from the symmetric limit represent single particle 
states. From Fig. 14 we can see that a linear extrapolation to of these two points would give consistent results with the 


GMO fit. The errors on the final values in Table |VII] are of similar sizes or larger than our estimate of the size of the discretization 
errors, in the range of 10 — 20 MeV. Apart from the positive parity singly charmed baryons, the shift AX^uHipiet does not have 
a significant impact on the results. 


Particle Parity 

Mi( 

SU(4) 

GeV) 

HQET 

M2 

SU(4) 

(GeV) 

HQET 

Sc H“ 

“C + 

Q.C + 

2.414(13)(22) 
2.534(09) (20) 
2.642(07) (18) 

2.434(20)(26) 

2.543(12)(21) 

2.648(09)(19) 

3.090(94)(95) 

3.177(67)(69) 

3.268(47)(50) 

3.144(115)(116) 
3.226(72)(74) 
3.294(50)(53) 

Ac + 

+ 

2.280(15)(23) 

2.442(11)(21) 

2.280(17)(24) 

2.442(11)(20) 

3.183(75)(77) 

3.234(37)(41) 

3.181(76)(78) 
3.245(36) (40) 

“cc 

^cc 

3.610(09)(12) 

3.713(06)(10) 

3.610(09)(12) 

3.713(06)(10) 

4.017(39)(40) 

4.157(31)(32) 

4.018(40)(41) 

4.158(31)(32) 

s* + 

H* + 

Uc -|- 


2.506(18)(25) 

2.608(13)(22) 

2.709(11)(21) 


3.234(121)(122) 
3.295(87) (89) 
3.355(64) (66) 

+ 

Ucc + 


3.694(07)(11) 

3.785(06)(10) 


4.078(36)(37) 

4.215(26)(28) 

E, 

flc 

2.698(87) (89) 
2.850(44) (47)“ 
2.995(20)(26)“ 




A, 

“c 


2.578(144)(145) 

2.761(77)(79) 



“cc 

^cc 


3.892(47)(48) 

4.061(28)(29) 

4.306(45)(44) 

4.394(42)(41) 


“c 


2.740(48)(51) 

2.891(32)(36) 

3.016(32)(37) 



'—'cc 


3.989(58)(58) 

4.132(42)(43) 


4.447(57) (57) 
4.567(55)(56) 


TABLE VII. Mass estimates at the physical point for the ground and first excited states of the singly and doubly charmed positive and negative 
parity baryons. The mean is computed using Eq. for the fits indicated in Table pg| and includes the shifts AAmuitipiet- The first error 
is due to varying the fit parameterization and is calculated using Eq. while the second error includes the uncertainty associated with 
AXniuitipiet. Note that for the negative parity channels we only give results for the operators for which the mass could reliably be extracted. 
The superscript ** indentifies likely scattering states. The results for the first excitations may also contain such states. 


Finally, we present the extrapolations for the mesons. The charmonium channels are fitted individually as for the other flavor 
singlets discussed above (to a constant and a constant plus quadratic term). Figure [T6| displays the extrapolations for the heavy- 
light mesons. Table |XII| (Appendix |E|l details the fit functions used. Since the heavy-light mesons form triplets, the expressions 
Eqs. ( [q] ) and ( [id] ) apply. While the lowest lying 0” and 1“ states exhibit GMO dependence on 5mi, this does not hold for the 
higher states. In particular, in the O'*' channel, the lowest eigenvalue for the c£ interpolator decreases rapidly towards 
and is very close in value to the corresponding non-interacting D-k state on each ensemble. In contrast, replacing the light quark 
with a strange quark, the Z/*g is below DK and is consistent with experiment at the physical point. Unfortunately, we were 
not able to reliably extract higher eigenvalues in these channels, nor did expanding the interpolator basis to include 0745 prove 
successful. Further work, including interpolators with derivatives will be performed in the future. A similar picture emerges for 
the 1+, however, here we are able to resolve closely lying levels by using two interpolators (see Table [n|. For the D meson case, 
the lowest eigenvalue coincides with D*tt, while the next level is much less sensitive to the light quark mass and is compatible 
with experiment at the physical point, leading us to associate this level with the Di. Switching to the Dg sector, the lowest 
level is below the D*K threshold and is consistent with the £)si(2460). The next level is both consistent with D*K and with 
Z/si(2536). A careful finite volume study is needed to clarify this issue. 
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FIG. 14. Extrapolation to for (left) positive parity and (right) negative parity charmed baryons on 1/ = 32^ x 64 ensembles. The 

final (extrapolated) results are displayed as in Fig.^ with the shift AXmuitipiet included. The type of fit displayed (GMO or LC) corresponds 
to that used to determine the final values, see Table ^ For improved visibility, where GM02, LCl and LC2 fits were used, we only show 
the GM02 fit. The type of interpolators chosen for each multiplet (HQET/SU(4)) is given above each plot. Where available the experimental 
results are displayed on the very left. 
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D. Mass differences 


Mass differences uncover specific aspects of the quark dynamics within the charmed baryons. We focus on three types 
of positive parity ground state mass differences, spin (J = |) — (J = |), “structure” (between the sextet and anti-triplet 
multiplets) and flavor (within multiplets corresponding to changing s u/d and between multiplets for c —> u/d and c —> s). 
The splittings are computed from the jackknifes of the individual masses and are extrapolated to using functional forms 

derived by taking the differences of the expressions for the masses. The parameterizations used for each channel are listed in 
Table IXIIII With the assumption that the singly charmed multiplets all require shifts, AXmuitipiet, of a similar size to make 
contact with the physical point, only the c ^ u/d and c —s flavor splittings require a correction. The final results are displayed 


in Fig. 17 


• Spin splittings: according to the predictions of HQET/pNRQCD, the | — 5 mass differences should be small, of 

0(A^/toc), vanishing in the static limit. Figurej^shows all spin splittings are indeed small, in the region of 60-85 MeV. 
They are also similar in size, independent of the light quark content and the number of charm quarks. This can be un¬ 
derstood as follows. For doubly charmed baryons, if the two charm quarks form a spatially small diquark (r < A ), 
the spin difference should correspond to 3/4 times the corresponding D (or Ds) fine structure splitting, where a factor 
3/2 is from coupling a = 1 or = 0 diquark, rather than a spin-1/2 antiquark, to the light quark. Another factor 
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FIG. 16. Extrapolations to the physical point for (top left) the D and Ds mesons, (top right) the Dq and D*q mesons, (bottom left) the lowest 
two eigenvalues for Di and (bottom right) the lowest two eigenvalues for Dai. For simplicity one fit is displayed in each case. The full set of 
fits used to calculate the final values are given in Table |Xli| (Appendix. These values include the shift X^uitipiet (and corresponding error) 
determined from the spin and flavor averaged D/Ds mass. For the Dq, Di and Dsi the results for the relevant non-interacting two particle 
states are also shown. 


1 /2 is due to the different color contraction within a baryon, relative to the meson color trace Il63lp^ Testing this HQET 
picture (or, equivalently, the short-distance pNRQCD picture), we have |(Md» — Me,) = 95(8) MeV compared to 
= 85(9) MeV and |(Md. - Md,) = 92(2) MeV compared to - Mn,, = 71(4) MeV. The Vt^c 

splitting agrees less well with this expectation, however, higher order corrections in HQET and (p)NRQCD of size A / 
and rricV^, respectively, are not necessarily small at the charm quark mass. Brown et al. Il43l . making a similar compar¬ 
ison (using the experimental meson splittings), find that agreement improves for bottom quarks. We remark that singly 
charmed baryons, where the two light quarks form a diquark, share the same spin- and color-factors of doubly charmed 
baryons (§•5 = 1)^ relative to charmed mesons. Therefore, also the Sc and flc fine structure splittings are of similar 

sizes, see Eig.[^ Finally, for r > A ^, doubly charmed baryons become similar to charmonia. Again, the spin and color 
factors are the same but in this case one should compare to the charmonium fine structure, |(Mj/^ —M^^) = 76(1) MeV. 
The situation is now reversed with the flee splitting being in closer agreement. Therefore, based on the fine structure, it is 
not possible to cleanly discriminate between the two cases (HQET or charmonium-like). As discussed in Section]^ from 
the meson sector we expect discretization effects in spin splittings of about 10-20 MeV, while the final errors in the figure 
are around 5-10 MeV. Nevertheless the sextet differences are reasonably consistent with experiment. 

• Flavor structure splittings: these are differences between baryons with the same and quark content but different 
flavor structure. The Ac (/ = 0) and Ec (/ = 1) lie in the anti-triplet and sextet multiplets, respectively, and differ in 


10 


In pNRQCD this factor is due to the light quark interacting with two color charges, in HQET one would ascribe it to the mass of an effective Q = cc 
antiquark. 
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FIG. 17. The spin and flavor mass splittings for the positive parity charmed baryons on the V = 32^ x 64 ensembles. For the splittings where 
the charm quark is replaced by a strange or light quark (bottom plot) two errors are shown indicating the uncertainty with and without the error 
on AXmuitipiet; this shift cancels in the other splittings. The experimental values are indicated, where we take the assignment suggested 
by the PDG in the cases where the quantum numbers have not been identified. The experimental errors are shown as bands. For those splittings 
for which a value is not quoted by the PDG, the error is calculated adding the uncertainties of the individual masses in quadrature. 


isospin. Similarly, the 2^ and the 2c share the same valence quark content but differ in terms of their SU(3) wavefunctions. 
Nature respects isospin symmetry reasonably well but not so SU(3); the physical 2^, and 2c states may be mixtures of the 
sextet and the anti-triplet states. While we have not studied this mixing we note that Brown et al. Il43l found that such 
effects are not significant. In Fig.[T7]one can see that the splittings are compatible with experiment. 


• Flavor splittings: in the middle panel of Fig. 17 we show differences between masses of particles where one or two 
strange quarks are replaced by light quarks. The mass differences for ss it are roughly twice as large as those for 
s ^ 1. Agreement is found with experiment for differences involving J = \ singly charmed baryons, however, our 
results seem to be systematically (slightly) lower for J = |. In the last panel of Fig. [lAwe show differences between 
masses of doubly charmed baryons and their singly charmed counterparts, replacing eitnCT one charm quark by a strange 
quark or by an up/down quark. Clearly, Mq. — Mq* « — Mn,, since « Mq. — etc., as already 

observed above under “Spin splittings". However, also and — M-bj^ « 

suggesting the dynamics of doubly charmed baryons to be closely linked to those of singly charmed baryons (and charmed 
mesons); in the first case the two heavy quarks seem to form a diquark (meaning r ^ A ) while in the second case either 
the light quarks form a diquark, interacting with the remaining charm quark or one light quark interacts with a heavy-light 
core. This would explain why the two systems resemble each other closely, up to constant shifts, due to differences of 
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the respective (constituent) quark masses. Note that replacing the remaining charm quark by a light quark flavor will not 
result in constant energy shifts. For example, the fine structure of the light and heavy baryons are very different. 


E. Comparison with previous results 

Our final results for the individual masses are summarized in Fig. where for the negative parity channels only the results 
for interpolators which gave a reliable signal are shown. Overall, for the singly charmed baryons there is reasonable agreement 
with experiment, as also seen for the mesons with open and hidden charm, which is encouraging in terms of the size of the 
remaining systematics. In particular, in some previous calculations Mcqq — or Mccq — splittings were computed, 

instead of the individual baryon masses, in order to remove the dependence on the charm quark mass and to reduce the size of 
discretization effects and uncertainties arising from tuning the charm quark mass using the meson sector. This has not been done 
here. However, states involving strange quarks are systematically lower than experiment, in particular, the He. This may be due 
to residual effects related to simulating along an incorrect fn = constant line. Furthermore, our finite volume study suggests 
that the positive parity ground state singly charmed baryons and excited state doubly charmed baryons will decrease by a few 
tens of MeV in the infinite volume limit. 
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FIG. 18. Singly charmed (top) and doubly charmed (bottom) low lying baryon spectra. (Left) our results (RQCD) from the V = 32^ x 64 ensembles. (Right) 
a comparison with previous determinations from Na et al. ED, Liu et al. Briceho et al. Ezi, Dlirr et al. do), ILGTI (39), PACS-CS collaboration gi), 
ETMC (42), Brown et al. (43) and HSC (B). Note that our results for the negative parity, spin-1/2 Qc and are likely scattering states, see Section [Tv| The 
first excitations may also contain such states. 


There have been a number of recent lattice studies of the charmed baryon spectra, also displayed in Fig.[T^ most of which are 
restricted to positive parity ground states. Different systematics apply in each case depending on, 

• the number of sea quarks included: Nf = 2 (Diirr et al. (401 1. Nf = 2-1-1 (Na et al. (TSl . Liu et al. PACS-CS iHdl . 
Brown et al. ED, HSC (Bl) and V/ = 2 -f 1 -f 1 (Briceno et al. lEll, ILGTI (ID, ETMC (421 1. 
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• The charm quark action; the Fermilab action Ea (Na et ah, Liu et al.), a relativistic heavy quark action GHEa (Briceno 
et al., PACS-CS, Brown et al.), twisted mass ll76l (ETMC), anisotropic clover IfTTlI (HSC), the Brillouin action HO) (Diirr 
et al.) and the overlap action (ILGTI). 

• The light valence and sea quark action: MILC Asqtad staggered sea 123 and valence (Na et al.), MILC Asqtad staggered 
sea ||79]| with domain wall valence (Liu et al.), MILC HISQ sea lISOl with tadpole-improved clover valence (Briceno et 
al.), QCDSF non-perturbative clover sea ISTl with Brillouin valence (Diirr et al.), MILC HISQ sea 1801 with overlap 
valence (ILGTI), non-perturbative clover sea l8^ and valence (PACS-CS), twisted mass sea l42ll and valence (ETMC), 
RBC/UKQCD domain wall sea l83]l and valence (Brown et al.) and anisotropic tadpole-improved clover sea and 
valence (HSC). 


Systematics arise from finite lattice spacings, unphysical quark masses, finite volumes, excited state contamination and, among 
mixed action approaches, violations of unitarity. Notable are the Briceno et al. El, ETMC na and Brown et al. El 
studies, involving both continuum and chiral limit extrapolations, the other works are predominately at a single lattice spac¬ 
ing. Despite such varied approaches, there is general agreement for the spectra. In particular, the lattice results are ap¬ 
proximately 80 MeV above the SELEX measurement of = 3518.7(1.7) IfTSl . The recent result of Borsanyi et al., 

M~++ — M~+ = 2.16(11)(17) MeV Il84l also contradicts the SELEX value of approximately 60 MeV for this isospin split- 
tinFIni. 

Of particular interest for our work are the determinations of ground state negative parity particles by ILGTI ll^ and the 
ground state and excitations for both parities for the doubly charmed spectrum by HSC Il44ll . ILGTI use pion masses down to 
240 MeV and two lattice spacings, a = 0.06 fm and 0.09 fm, both with a spatial extent L = 2.9 fm while HSC employ a 
temporal lattice spacing of at ~ 0.035 fm on anisotropic lattices with L = 1.9 fm and a single pion mass, — 391 MeV. 
Considering the different systematics involved in each study, Eig.[T8| shows that the determinations for the ground state negative 
parity channels are reasonably consistent, including those which we identify as scattering states. The picture is less satisfying 
comparing with the HSC predictions for excited doubly charmed baryons. Our positive parity first excited states are significantly 
below those of HScj^ even for the . This may be due to finite volume effects and/or the lack of chiral extrap olation in 


the HSC study. Certainly the discrepancy is smaller for the compared to the See’’ and, as discussed in Section 
expect finite volume effects to increase the mass for smaller physical volumes. Eor the negative parity states, HSC predict several 
levels (as expected from non-relativistic quark models with SU(6)xO(3) symmetry, see Ref. Il44l for a discussion) in the range 
4.0 — 4.8 GeV. With only one interpolator (but multiple smearings) we are not able to resolve closely lying states and our excited 
states are above the first three HSC levels. With our limited basis the first excitation could also have an overlap with a state 
corresponding to higher continuum spin. This would correspond to J = | for spin-1/2 interpolators and J = | for spin-3/2. A 
more extensive analysis is needed in order to make a closer comparison. 


^(*) 


VIA 


we 


VII. CONCLUSIONS 

In this work we have presented results for the ground states and first excitations of singly and doubly charmed baryons. 
Through the use of the variational method with a basis of three differently smeared interpolators, reliable signals were obtained 
for both positive and negative parities, where for the latter the first excitation was only extracted for doubly charmed states. Eor 
spin-1/2 channels we implemented both HQET and SU(4) inspired interpolators and found consistent results for positive parities, 
while for negative parities SU(4) interpolators had a much better overlap with the ground state for sextet baryons (Ec, 5),, flc) and 
HQET interpolators for the anti-triplet (Ac, Sc) and triplet (Sec, ^cc)- Overall agreement was obtained with experiment, where 
available, suggesting the remaining systematics are not large. This is supported by the results for the lower lying charmonium, D 
and Dg mesons, which also reproduced the gross experimental spectra. Discretization effects can be significant in a simulation 
at a single lattice spacing, for which arric ~ 0.4. Eine structure is particularly sensitive to such effects, and in the meson sector 
spin splittings were underestimated by 10-20 MeV. However, for the singly charmed baryons the spin splittings are consistent 
with experiment within statistical uncertainties ranging from 5 to 11 MeV. Our study of finite volume effects suggests the lightest 
members of the ground state singly charmed multiplets and first excitations of the doubly charmed triplet, in the infinite volume 
limit, could be around Icr lower (in terms of the final quoted error). 

The negative parity channels and first excitations of both parities require careful study due to a number of decay channels with 
thresholds close to our mass estimates. By looking for degeneracy of the masses with the threshold values on each ensemble, 
over the range of pion masses = 259 — 461 MeV, we identified the ground states in the flc and S), negative parity channels 
as scattering states. Interestingly, for the singly charmed anti-triplet the HQET interpolators had a good overlap with the level 
consistent with a bound state, while the SU(4) interpolators gave results consistent with the corresponding thresholds. We 


*' For improved visibility we have omitted the HSC results for positive parity higher excitations from Fig. 18 
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also encountered scattering states in the meson spectra, where, by including additional interpolators in the variational basis, 
results consistent with experiment for the Di and Dsi were obtained. Cross-correlation functions between the SU(4) and HQET 
interpolators would be needed to repeat this analysis for the baryons. 

One of the main aims of our study was to investigate the light flavor structure of charmed baryons; employing ensembles 
with the average light quark mass fixed enables expansions in the flavor symmetry breaking quark mass difference, starting 
from a point where the lattice results are most precise. The cotTesponding Gell-Mann-Okubo relations worked well for almost 
all multiplets and the extrapolations were well constrained, with the linear terms in 5mt within each multiplet being fixed by 
one coefficient. In some cases the quadratic terms were statistically significant, however, the magnitudes of these contributions 
were still small. A notable exception was the singly charmed anti-triplet (Ac, Sc), which did not fit the GMO pattern. This 
may indicate a more complex internal structure. However, one can still perform a phenomenological fit using an expansion in 
&mi, where the linear coefficients are not constrained. The errors on our final values were estimated conservatively, including 
both linear and quadratic fits. An additional uncertainty is present in our results due to the simulation trajectory missing the 
physical point in the mass plane. This effect was corrected for in a consistent way and the agreement found with experiment 
for extrapolated flavor singlet mass combinations suggests the correction is small. However, we note that our results for heavy 
baryons involving strange quarks are consistently lower than experiment which may indicate residual effects. 

The dynamics of quarks within heavy baryons is interesting in terms of possible similarities with the heavy-heavy or heavy- 
light mesons and the applicability of various effective theories. In particular, for doubly charmed baryons, there are two relevant 
limits, in terms of a point-like heavy-heavy diquark (the HQET or short distance pNRQCD picture) or a more diffuse heavy-light 
diquark (charmonium picture). Elavor splittings where the charm quark is replaced by a strange or up/down quark suggest the 
dynamics of doubly charmed baryons are closely linked to those of singly charmed ones. 

In the future we plan to perform a more extensive study on CLS configurations ifSOll . with a larger range of pion masses, all 
with LM,r ^ 4. This will enable a further investigation of the magnitude of flavor symmetry breaking effects. Initial studies 
show the physical trajectory has been reproduced reasonably well ll85l . 
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Appendix A: Meson effective masses 


Effective masses of the lowest two eigenvalues for a number of channels for charmonium and D/Dg are shown in Fig. 19 for 
the larger symmetric ensemble. The effective masses for the lowest four eigenvalues for Di and Dgi on the = 259 MeV 


ensemble are displayed in Fig. 20 
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FIG. 19. Sample of effective masses of the charmonium and D/Ds mesons for the symmetric ensemble with V = 32^ x 64. The filled 
regions indicate the fitting ranges chosen and the fit result including errors. 
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FIG. 20. As in Fig. 19 for the effective masses of the lowest four eigenvalues (EV) for the 1’*' channels on the = 259 MeV ensemble 
with V = 32^ X 64. 
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Appendix B: Finite volume effects 


The finite volume effects computed using the matched V = 24^ x 48 and 32^ x 64 ensembles are given in Tables 


Vlll 


and 


IX for positive and negative parity particles, respectively. 


Particle 

p 



Eirst eigenvalue 



Second eigenvalue 





SU(4) 

HQET 

SU(4) 

HQET 




AM (MeV) 

ZiAi //W X 

AM (MeV) 

A Ai /c/ \ 

AM (MeV) 

AfPf (OJ \ 

mP, 

AM (MeV) 


Ec, H),, Qc 

+ 

f^sym 

33(10) 

1.3 

24(11) 

0.9 

-63(89) 

-1.9 

-59(81) 

-1.8 

Ec 

+ 


31(15) 

1.2 

15(16) 

0.6 

69(134) 

2.1 

-4(139) 

-0.1 

' —'C 

+ 

^asym 

18(11) 

0.7 

16(11) 

0.6 

50(102) 

1.5 

-19(96) 

-0.6 

Qc 

+ 


10 (8) 

0.4 

8(9) 

0.3 

31(84) 

1.0 

-22(77) 

-0.7 


+ 

f^sym 


39(12) 

1.5 


-139(90) 

-4.0 

E* 

+ 




27(17) 

1.0 



129(164) 

3.8 

'—'C 

+ 

^asym 



15(12) 

0.6 



67(108) 

2.0 

q: 

+ 




6 (10) 

0.2 



34(84) 

1.0 

Ac, He 

+ 

f^sym 

32(8) 

1.3 

33(8) 

1.4 

-229(98) 

-7.0 

-266(98) 

-8.1 

Ac 

+ 


19(12) 

0.8 

20(13) 

0.9 

-139(187) 

-4.4 

-150(201) 

-4.7 


+ 

^asym 

11(9) 

0.5 

12(9) 

0.5 

-63(121) 

-1.9 

-72(126) 

-2.2 

“CC 

+ 

^sym 

4(6) 

0.1 

3(6) 

0.1 

76(24) 

1.9 

79(24) 

1.9 

' —'CC 

+ 


4(8) 

0.1 

6(7) 

0.2 

63(29) 

1.5 

56(29) 

1.4 

^CC 

+ 

^asym 

-1(6) 

-0.0 

0(5) 

0.0 

43(27) 

1.0 

39(27) 

0.9 

O* 

^ “CC 

+ 

f^sym 


9(7) 

0.3 


87(26) 

2.1 

'—'CC 

+ 




1(9) 

0.0 



69(30) 

1.7 


+ 

^asym 



-5(7) 

-0.1 



46(27) 

1.1 


TABLE VIII. Mass differences, AM — M 24 — M 32 , between the positive parity masses extracted from the V = 24® x 48 and the 
32® X 64 ensembles, where the quark mass parameters (k) are the same: the symmetric point Ksym = 0.1209 and the asymmetric combination 
Kasym = (0.12104, 0.12062). The errors on the differences are computed by adding the individual uncertainties in quadrature. 


Particle P 

Operator 

AM (MeV) (%) 

First eigenvalue 

^Ci “c? ^c “ ^sym. 

SU(4) 

-36(37) -1.2 

Ec 

“c “ ^asym 

He 

SU(4) 

-49(40) -1.7 

-37(31) -1.3 

-22(24) -0.7 

A.C, ~ ^sym 

HQET 

-4(40) -0.1 

Ac 

/^asym 

“C 

HQET 

122(48) 4.7 

63(36) 2.3 

i^cc? “CC “ ^sym. 

HQET 

9(23) 0.2 

“CC 

Q /^asym 

“ “CC 

HQET 

-2(23) -0.1 

-8(17) -0.2 

^cc? “CC " ^sym. 

HQET 

-19(23) -0.5 

—CC " 

f^asym 

HQET 

-39(22) -1.0 

-31(16) -0.8 


Second eigenvalue 


l^cci “CC “ f^sym 

SU(4) 

22(47) 0.5 

“CC 

Q ^asym 

i “CC 

SU(4) 

43(46) 1.0 

32(39) 0.7 

Q* '=* _ ^ 

^ ^CCI “CC ^sym 

HQET 

-67(39) -1.5 

“CC " 

'^asym 

HQET 

79(52) 1.8 

21(36) 0.5 


TABLE IX. As in Table |VIlI| for the negative parity states. 
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Tensor 

rank 

Representation 

S 

(0,0) 

1 

q- 

(1,0) 

3 

Qa 

(0,1) 

3 

ot 

(1,1) 

8 

gab 

(2,0) 

6 

Sab 

(0,2) 

6 

j-^abc 

(3,0) 

10 

Dabc 

(0,3) 

10 

rpcd 
^ ab 

(2,2) 

27 


TABLE X. Tensors corresponding to the lower dimensional irreducible representations of SU(3). 


Appendix C: Tensorial notation for SU(3) representations 


For convenience we review the tensor notation for SU(3) representations. Consider a general SU(3) tensor with n 

contravariant indices and m covariant indices (i.e. rank (n, m)), which transforms according to 


rpdx.. .CL 
^bi...byr 




(Cl) 


where U is the SU(3) transformation. We can construct lower rank tensors starting from by performing the following 

contractions. 


cbj >T^CLl . . .dn 

^ai^bi...bm. ’ 

for 

1 < 2 < n, 1 < j < TO, rank (n — 1, to — 1), 


rriClX .. .dn 

iaiajbm+l^bi...br,^ 

for 

rank (n — 2, TO + 1), 

(C2) 

bibjdn + 1 'Y'Q-l ■ ■■O.n 

^ ^bi...bm 

for 

1 < *,J < w, rank (n + 1, TO — 2), 



where 6^, Cabc and are SU(3) invariant. Irreducible tensors cannot be expressed in terms of tensors of lower rank, and hence, 
one obtains zero from any of the above contractions. This means that the irreducible tensors are pairwise symmetric in both the 
covariant and contravariant indices and traceless on any index. The lower dimensional irreducible tensors are listed in Table [X] 
with the dimension D{m, n) = + l)(m + l)(n + to + 2). We are interested in decomposing the product of two irreducible 

representations. The tensor product is expressed in terms of the irreducible tensors and 6^, tabc or For example, 

303 = 1 ©8, q'^qb=^5tS + Ol. (C3) 

Similarly, tensors of irreducible representations can be expressed in terms of the components of the fundamental representation, 

S = q‘'qc and = q°^qb - (C4) 


Appendix D: Gell-Mann-Okubo relations for charmed baryons 

Below, we derive the expressions for the mass dependence of each state within a multiplet up to hrst order in the flavor 
symmetry violating parameter 5mi. In Eq. (0 we have matrix elements of tensor operators between tensor states 

{Br\H\Br), (D1) 

where the baryons B and the Hamiltonian H are in the irreducible representations R and Tr, respectively. Of interest are 
i? = 6, 3, 3 and Tr = 1, 8. All the indices of the representations must be contracted, i.e., in combination they form the 
trivial representation. Thus, the number of terms each matrix element can depend on is equal to the number of times the singlet 
representation appears in the decomposition of the direct product R 0 Th 0 For the combinations of interest, illustrated 
below, there is only a single instance of 1 in the direct products. 

In the following, for each representation we first consider R 0 R and then the direct product with Tr. 


12 


This is a generalization of the Wigner-Eckart theorem. 
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1 ± Q± 

• Sextet; corresponding to the singly charmed baryons with J — ^ or J = | . We have 

R(g)R = 6(g)6 = l©80 27 (D2) 

The 6 representation is constructed from components of the fundamental representation, 5“^ = (q°‘q^ + where 

= (u, d, s). Taking the example of the J = \ sextet, we can assign the singly charmed baryons to the elements of 

gab. 


511 = 

2uu - 2E++, 


522 

= 2dd 

^2E0, 

(D3) 


= {ud + du) ^ 

^V2X, 


= 532 

= {ds + sd) ~ 


^13 = 

- = {us + su) ' 

- V 2 X, 

533 

= 2ss ‘ 

-20°. 



The tensor product of 5“^ 
and 27-plet (T^^). 

6(g)6 : 

with 

S = SefSf^ 

Ol = SaeS^^ - (D5) 

= SabS^'^ - i {<5^(5,+ SiiSbeSn + S^,{SaeS^<^) + S^iSaeSn} + ^ i^a^b + S, 


with its conjugate representation, Scd can be decomposed in terms of a singlet (S), octet (O^) 


SabS^‘^ = 4 {S^st + SX} S + i {6Xt + sXt + 5X1 + 5toi) + t; 


cd 
ab 5 


(D4) 


where the summation over repeated indices is understood. 

We can now identify the singlet representation that appears in the product 6 i 
elements. For Th = 1 we simply obtain the singlet term appearing in Eq. (D4i. Thus, the lowest order term in the 
expansion in 5mi (cf. Eq. Q) is given by 


Th © 6 and the expressions for the matrix 


m{B-^\Ho\B^) = |E++|2 + |E+| 


^0|2 


l'='^+|2 _i_ |'^^0|2 


io: 


0|2 


(D6) 


The coefficient is the mass of the sextet in the flavor symmetric limit. 

Eor the next order term, iFg cx Tg cx Of, where Tg is the second diagonal SU(3) generator and O is the octet represen¬ 
tation (Th = 8 = 8). In the direct product 60806 = 80(1 ©8© 27) the singlet appears in the 8 0 8 term. More 
explicitly, one contracts Og from the Hamiltonian with the = S^cS'^^ — XcdS‘^‘^ part of the tensor product of the 


baryon representations in Eq. (D4l; 


Sme{B-e\Hs\Be) = [-| (|S++|2 + jE+p + XX) + | (| 


7T -t|2 


|~' 0|2 




(D7) 


Note that all particles in the sextet have the same coefficient, Ai, for this 0{5mg) term modulo factors that are proportional 
to the light hypercharge. Isolating the mass terms for each particle up to first order in the Taylor expansion and taking the 
isospin limit (i.e. dropping the superscripts for the electric charges), we arrive at 


= Mq - + 0{5mg), 

Ms, = Mg® + + 0(,5 to^), (D 8 ) 

Mq, = M® + Xi5mg + 0{Smj). 


The same expressions, but with different coefficients in each case, can be used for the other sextets. 

• Triplet; corresponding to the doubly charmed baryons with J = ^ or J = | . We proceed in an analogous way to the 
above starting with R0R = 303 = 8©1. The assignment of the components of the fundamental (3) representation 
to the corresponding doubly charmed baryons is straightforward. Eor the example of J = 4 , 

= S ^ H+. 


q^ = u 


"-I-+ 

“CC f 


q^ = d 


"+ 

^CO 


(D9) 
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The tensor product of ( 7 “ with the conjugate representation qh takes the form 

3® 3: q^qb = ^6^S + Ol (DIO) 

For the lowest order matrix element, corresponding to Th = 1, we find 


m{B-s\Ho\Bs) = [| 


3 r I'1:7++12 I 177+12 


\^tc? 


corresponding to the singlet appearing in Eq. (DIO 1 , while at the next order with Th = 8 

6mi{Bs\Hsm = 5m,B^ [-|(|S++p + |S+H + ||f]+p] ^ 


(Dll) 


(D12) 


arising from contracting 0\ from the Hamiltonian with (9| = q^q-^ — ^q'^qc from Eq. (DIO 1 . 

Considering the individual states we arrive at the following Gell-Mann-Okubo relations for the doubly charmed baryons, 


= Mq - ^BiSm, + O(Smj), 

= Ml + lB,5m, + 0{6ml). (D13) 


As before, the same expressions apply to the other triplets but with different coefficients. 

• Anti-triplet; cori'esponding to the singly charmed baryons with J = . Here R 0 R is the same as for the triplet. We 

construct the 3 representation from the 3, qc = ecab<l°‘q^, leading to the components, 

qi = {ds — sd) ~ 5°, q 2 = {us — su) ~ q^ = {ud — du) ~ A+, (D14) 


in terms of the corresponding singly charmed baryons. Using this and the tensor product, q°'qb, given in Eqs. (DlOi we 
can derive the expressions, as above, for the Taylor expansion 


= Mg — ^CiSm, + 0((5m^), 

Mb., = Ml + lCi6m, + 0(<5m^). (D15) 
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Appendix E: Fitting and extrapolation details 


The fi tting ranges chosen to extract the charmed baryon masses and the charmonium, D and Ds spectra are given in Tables |XI 
and 


XII 


respectively. Also included are the fits employed in the extrapolation to that are used to compute the physical 


point result, see Eqs. (fT7]i and (H^. Similarly, for the baryon mass differences in Table IXIIII In Table |XIV1 we provide the 

1io( 


coefficients extracted from linear and quadratic fits using the GMO formulae (Eqs. m to (|l6|l) and their less constrained 


counter-parts (denoted LCl and LC2) for the multiplets appearing in Eig. 14 


Particle P Operator 

First eigenvalue 

Fit-range Extrapolation 

Second eigenvalue 
Fit-range Extrapolation 

Sc + 

S'c + SU(4) 
r^c + 

[10 — 24]Nf=3, [10 — 18|M,r=255 MeV 

24] GM02, LCLLC2 

[9-13] GM01,GM02 

+ SU(4) 

—c H” 

[11-22] GMOl, GM02 

[7-13] GMOl, GM02 

t SU(4) 

i ^CC 1 

[11-24] GM02, LC1,LC2 

[10-20] GM01,GM02 

Sc + 

+ HQET 

r^c + 

[10 - 19] GMOl, GM02 

[9-16] GMOl, GM02 

+ HQET 

— C H” 

[11-22] GM01,GM02 

[7-13] GMOl, GM02 

““ + HQET 

i ^CC 1 

[11-25] GM02, LCLLC2 

[10-20] GM01,GM02 

s* + 

H* -h HQET 
H* -I- 

[10 - 20] GMOl, GM02 

[9-13] GMOl, GM02 

+ HQET 

^ ^cc 

[10 - 25] GMOl, GM02 

[10-20] GM01,GM02 

Sc - 

H'c - SU(4) 

He - 

[10 - 15] GMOl, GM02 


Hec - SU(4) 


[8-13] GMOl, GM02 

~ HQET 

“C 

[11-15] LC1,LC2 


“ HQET 

i ^CC 

[10-16] GM02,LC1,LC2 


s; - 

H* - HQET 
He - 

[10 - 16] GMOl, GM02 


“ HQET 

^ ^CC 

[9-15] LC1,LC2 

[7 - 13] LCl, LC2 


TABLE XL Fitting range s used to extract the charmed baryon masses and the fit functions employed for computing the final results at the 
physical point (see Section VIC ' for the 32^ x 64 ensembles. Note that for the negative parity states we only give details for the operators for 


which reliable signals could be obtained. 
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Particle 

First eigenvalue 
Fit-range Extrapolation 

Second eigenvalue 
Fit-range Extrapolation 

^c cr^ 

j/^ 1— 

XcO 0++ 

Xcl 1 + + 

he 1+- 

15-23 FSLFS2 

15-23 FSLFS2 

15-27 FS1,FS2 

15-27 FSLFS2 

13-27 FSLFS2 

10-17 FS1,FS2 

10-17 FS1,FS2 

O 

1 

15-26 GM02, LCl, LC2 


D* 

(d:) 

15-26 GMOl, GM02 


a Q + 

cs (Dto) ^ 

NCI NC2 
10-15 ^ ^ 

NCLNC2 


(C7i75s), (ceijkXjJks) 

NC1,NC2 
NCI, NC2 

NC1,NC2 

NC1,NC2 


TABLE XII. As in Table |Xl| for the lower lying meson channels. The extrapolations ESI and FS2 refer to flavour singlet fits which include a 
constant term and a constant plus a quadratic (5mj) term, respectively. Similarly, NCI and NC2 denote a linear fit to a single channel and a 
linear plus quadratic fit, respectively. The GMO and LC fits are defined in the Section[VIB| 


Difference 

Operator 

Extrapolation 

Spin splittings 

y* y ■=* o* o 

HQET 

GMOl, GM02 

^CC “CC? ^cc ^cc 

HQET 

GM02, LCl, LC2 

Structure splittings 

V A 

—‘C “C 

SU(4) 

GMOl, GM02 

Ec - AcE'e - Ec 

HQET 

GM01,GM02 

Flavour splittings s —> u/d 

E'c - Ec, Oc - H'c, Oc - Ec 

SU(4) 

GMOl, GM02 

H'c - Ec, Oc - E'c, Oc - Ec 

HQET 

GMOl, GM02 

He - Ac 

SU(4) 

LC1,LC2 

He - Ac 

HQET 

LC1,LC2 

“cc 

SU(4) 

LC1,LC2 

“cc 

HQET 

LC1,LC2 

Q* Q* y* 

HQET 

GM01,GM02 

Flavour splittings c ^ u/d and c — ^ 

S 

O O - VO 

ilec Jtc ?'—CC “C ?'—CC ^Ci •^^CC “C 

SU(4) 

GMOl, GM02 

o O ~ TO =■' 

ii-cc 'itc?' —'CC “C ?'—CC “c 

HQET 

GM01,GM02 

O* o* V* O* 

—'CC “C ?'—CC ^C'l'^^CC “C 

HQET 

GM01,GM02 


TABLE XIII. The fit functions used for the extrapolation of the positive parity ground state differences to The splittings listed 

together were fitted simultaneously. Note that for He — Ac and Qcc — Hec, there are no corresponding GMO expressions, since these are 
splittings between multiplets containing two states. 
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Sextet SU(4) J = | + , (E^, H', n,) 


-^OtLCl) 

2.5530(59) 


-^0(LC2) 

2.5502(66) 

7l-l{LCl) 

-0.1428(153) 


j4i(lC2) 

-0.0921(601) 

^1{LC1) 

0.0232(127) 


^1(LC2) 

0.0720(497) 

^1{LC1) 

0.1740(115) 

LC2 

^1(LC2 

0.2155(422) 




^2(LC2) 

-0.0744(941) 




^3(LC2) 

-0.0715(723) 




^4(LC2) 

-0.0592(576) 


LCl 


GM02 


K(gmo 2) 2.5509(57) 
^i(GM02) 0.3006(282) 

^2(GM02) —0.0688(344) 
^3(GM02) —0.0382(181) 
^4(GM02) —0.0372(273) 


Sextet SU(4) = | , (E^, Q,) 


GMOl 

<(GMOi) 2.8650(160) 


'^0(GMO2) 

2.9080(239) 

j 4 i(gmoi) 0.3383(300) 


^1(GM02) 

0.0312(1065) 



GM02 

^2(GM02) 

-0.4212(1391) 




^3(GM02) 

-0.0875(734) 




^4(GM02) 

0.2048(1019) 



Triplet HQET = 

i"*" (=■ n i 

2 » 5 ^cc J 


LCl 

-^0(LC1) 

Sl{LCl) 

^1{LC1) 

3.6585(44) 

0.0723(89) 

-0.0729(99) 

LC2 

]\yf3 

^“0(LC2) 

^1(LC2) 

^1(LC2) 

^2(LC2) 

-B3(LC2) 

3.6596(48) 

0.0557(305) 

-0.0955(369) 

0.0242(420) 

0.0349(551) 

GM02 

'''^0(GMO2) 

5i(GM02) 

52(GM02) 

^3(GM02) 

3.6559(40) 

-0.0471(46) 

-0.0271(164) 

-0.0330(199) 








Triplet HQET = 

2 ’ (‘^CC7 ^cc) 



-^□(LCl) 

4.0143(191) 


-^0(LC2) 

4.0281(216) 


'^0(GMO2) 

4.0158(172) 

LCl 

Sl{LCl) 

0.0559(392) 


5i(LC2) 

-0.0817(1274) 

GM02 

5i(GM02) 

-0.0152(227) 


^1{LC1) 

-0.1580(445) 

LC2 

^1(LC2) 

-0.1611(1568) 


52(GM02) 

0.0435(759) 





52(LC2) 

0.1897(1730) 


^3(GM02) 

-0.2641(908) 





-B3(LC2) 

-0.0605(2348) 





Antitriplet SU(4) (Ac, He) 





GMOl 

-^OCGMOl) 

2.4274(40) 


'^0(GMO2) 

2.4337(56) 





Cl(GMOl) 

0.0756(25) 

GM02 

C'l{GM02) 

0.0817(90) 







C'2{GM02) 

-0.0099(537) 








C'3{GM02) 

-0.0429(205) 





Antitriplet HQET | 

[I] 

< 

1 







-^0(LC1) 

2.7825(332) 


-^0(LC2) 

2.8002(356) 

LCl 

CltLCl) 

-0.3500(936) 


Cl(LC2) 

-0.9128(3261) 


Cl(LCl) 

-0.0276(783) 

LC2 

Ci(LC2) 

-0.3978(2514) 




C2(LC2) 

1.0029(5636) 




C3(LC2) 

0.5979(3759) 


TABLE XIV. The coefficients determined from linear and quadratic GMO and LC fits (see Section [VIB^ to the ground state positive and 
negative parity multiplets displayed in Fig.|14| 
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